Dimensional crossover in quantum critical metallic magnets 
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Nearly magnetic metals often have layered lattice structures, consisting of coupled planes. In such 
a situation, physical properties will display, upon decreasing temperature or energy, a dimensional 
crossover from two-dimensional (2d) to three-dimensional (3d) behavior, which is particularly inter- 
esting near quantum criticality. Here we study this crossover in thermodynamics using a suitably 
generalized Landau-Ginzburg- Wilson approach to the critical behavior, combined with renormal- 
ization group techniques. We focus on two experimentally relevant cases: the crossover from a 2d 
to a 3d antiferromagnet, and the crossover from a 2d ferromagnet to a 3d antiferromagnet. As 
naive scaling does not apply at and above the upper critical dimension, two crossover scales arise 
which can be associated with separate dimensional crossovers of classical and quantum fluctuations, 
respectively. In particular, we find an intermediate regime with novel power laws where the quan- 
tum fluctuations still have a 2d and the classical fluctuations already have a 3d character. For the 
ferromagnet-to-antiferromagnet crossover, the mismatch of the dynamical exponents between the 
2d and 3d regimes leads to an even richer crossover structure, with an interesting 2d non-critical 
regime sandwiched between two critical regimes. For all cases, we find that thermal expansion and 
compressibility are particularly sensitive probes of the dimensional crossover. Finally, we relate our 
results to experiments on the quantum critical heavy-fermion metals CeCue-^Aurc, YbRli2Si2, and 
CeCoIng. 

PACS numbers: 71.27.+a,71.10.Hf 



I. INTRODUCTION 



Quantum phase transitions (QPT) in metals are a 
fascinating field of today's condensed matter research. 1 
Heavy-fermion materials play a prominent role: Fre- 
quently, one observes non-Fermi liquid behavior which 
is thought to be associated with an antiferromagnetic 
(AFM) instability of the itinerant electrons. The crit- 
ical spin fluctuations near the phase transition lead to 
unconventional power laws in transport and thermody- 
namic quantities at low temperatures. On the theory 
side, a Landau-Ginzburg- Wilson (LGW) description of 
the critical magnetic degrees of freedom, developed by 
Hertz; 2 - Moriya^, and Millis^ accounts for many of the 
experimental signatures of magnetic criticality. 

However, some heavy-fermion compounds do not eas- 
ily fit in the LGW picture in d — 3 spatial dimensions. 
For CeCu.6-xAu x , it was realized^ that most thermody- 
namic signatures of the QPT at x = 0.1 are consistent 
with the assumption of the underlying AFM spin fluctu- 
ations to be 2d. It came as a surprise that 2d spin fluc- 
tuations should prevail in an intrinsically 3d alloy, but 
the 2d character was subsequently confirmed in neutron- 
scattering experiments. 6 Importantly, the AFM order ob- 
served in CeCue-zAuz below the Neel temperature is 
fully 3d. On general grounds, one thus expects a di- 
mensional crossover within the paramagnetic phase from 
2d magnetic fluctuations at elevated temperatures to 3d 
fluctuations at lowest temperatures or in the immediate 
vicinity of the phase transition. Experimentally, the di- 
mensional crossover in CeCug-aiAua has proven to be 
elusive so far. 



A related heavy-fermion metal, not easily fitting the 
LGW theory framework, is YbRh2Si2^ It shows a phase 
transition to an ordered phase at 70 mK, which is be- 
lieved to be AFM, however, a confirmation by neutron 
scattering is not available to date. An additional as- 
pect is that YbRb.2Si2 seems to be almost ferromagnetic 
(FM), 8 and we will return to this issue later in this pa- 
per. The unusual properties of both CeCug-xAuz and 
YbRli2Si2 have prompted speculations on the inapplica- 
bility of the LGW theory, which describes a magnetic 
instability of well-defined quasiparticles: Instead, it was 
proposed that the Kondo effect, being responsible for the 
formation of the heavy quasiparticles, breaks down at 
the quantum critical point (QCP)^^ Different scenar- 
ios and theoretical descriptions of this Kondo breakdown 
have been put forward j 9 ' 10 i 11 The scenario of so-called 
"local quantum criticality"— uses an extension of dynam- 
ical mean-field theory to map the Kondo-lattice problem 
to a self-consistent impurity model where the Kondo ef- 
fect may be suppressed by critical bulk spin fluctuations. 
This particular scenario for a Kondo-breakdown QCP re- 
quires the spin fluctuations to be 2d; for 3d spin fluctua- 
tions this model predicts a conventional magnetic QCP of 
LGW type. As the critical spin fluctuations of the mate- 
rial can again be expected to become 3d at low energies, 
the local quantum criticality should be restricted to el- 
evated temperatures or energies above the dimensional 
crossover. 

Other examples of layered metals with magnetic QCP 
are the heavy fermions CeMIns (M=Co,Rh,Ir), the high- 
temperature superconducting cuprates and iron pnic- 
tides, and the metamagnetic ruthenate S^RusOy. 

The purpose of this paper is to study theoretically the 
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dimensional crossover of critical magnetic fluctuations in 
the framework of the LGW model, with focus on the 
2d AFM to 3d AFM and 2d FM to 3d AFM crossovers. 
Microscopically, we imagine a system of planes of inter- 
acting electrons, with tendency toward antiferromagnetic 
or ferromagnetic in-plane ordering, and a weak antiferro- 
magnetic inter-plane coupling. We will study how the di- 
mensional crossover is reflected in the correlation length, 
specific heat, thermal expansion, compressibility, and the 
Griineisen parameter^ Our primary goal is to identify 
observables which are suited for an experimental search 
for a dimensional crossover. 

In this paper, we shall refrain from a detailed mi- 
croscopic modeling of the magnetic inter-plane coupling. 
The underlying lattice geometry and bandstructure will 
influence some of the non-universal properties of the 
dimensional crossover; however, the existence of well- 
defined 2d and 3d regimes is independent of those details. 
We shall ignore complications arising from the possibility 
of the inter-plane coupling being frustrated ; 13 ! 14 ' 15 Even 
for a fully frustrated coupling, a dimensional crossover to 
3d behavior at low energies will generically occur, albeit 
with a possibly small crossover scale d 4 ' 15 (In situations 
with frustration, the effective 3d coupling within the or- 
dered phase may be enhanced due to order-from-disorder 
mechanisms.) 

We shall restrict our analysis to the framework of the 
LGW theory of itinerant spin fluctuations. It has been 
discussed that this approach may break down at low- 
est energies due to the occurrence of singular terms in 
the LGW expansion (i) for ferromagnets in both 2d and 
3dji£ and (ii) for 2d antiferromagnets^ These complica- 
tions will be ignored for simplicity, a justification being 
that the low-energy behavior of our model is invariably 
3d AFM (where the LGW approach is believed to be 
valid).— Our theory may also be combined with the ideas 
of local quantum criticality, 9 this is beyond the scope of 
this paper. Similarly, an explicit treatment of the mag- 
netically ordered phases shall not be performed herei 19 ' 20 



A. Summary of results 

The main results of our analysis are the crossover 
phase diagrams, Figs. [1] and H for the 2d - 3d AFM 
and the 2d FM - 3d AFM crossover, respectively. The 
large anisotropy in the spin-fluctuation spectrum defines 
a small momentum scale A,, [see Eq. below] that de- 
termines the positions of the dimensional crossovers in 
the phase diagram, as indicated by the shaded areas. 
Generally, the behavior changes from 2d to 3d upon ap- 
proaching criticality. 

However, the fact that the QPT under consideration 
are at or above their upper critical dimension renders 
naive scaling invalid; as a result, the quantum critical 
regime is characterized by two distinct length scales: The 
physical correlation length £ and a thermal length given 
by T~ 1//z , where T is the temperature and z the dynam- 




FIG. 1: Phase diagram of the anisotropic LGW model in the 
temperature-control parameter plane for the 2d AFM to 3d 
AFM crossover. The crossover scale is determined by the mo- 
mentum scale A,,. The 3d spin fluctuations dominate close to 
the quantum critical point. A dimensional crossover occurs 
upon increasing the distance to the QCP indicated by the 
shaded region. The phase boundary T c (r) changes its behav- 
ior at this crossover, see Fig. U below. There is an additional 
dimensional crossover at the temperature scale T c i(r) where 
£ ~ l/A v , associated with the classical critical fluctuations. 
The thin line, Tg(t), close to the critical temperature, T c (r), 
indicates the Ginzburg temperature where the crossover to 
classical Wilson- Fisher behavior occurs. The two lines T c i (r) 
and Tq(t) cross at a temperature T x . The dashed line sep- 
arates the low-temperature magnetically disordered (Fermi- 
liquid, FL) regime, T <C r, from the quantum critical regime, 
T > r. Note that the labels "2d" and "3d" refer to the behav- 
ior of the critical or near-critical spin fluctuations, for details 
see text. 




FIG. 2: Phase diagram of the anisotropic LGW model as in 
Fig. [TJ but for the 2d FM to 3d AFM crossover. Here, the 
critical 2d regime, T > or r > A^, is separated from the 
critical 3d regime, T < At, and r < A*, by a region where 2d 
non-critical Fermi liquid behavior prevails. 
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ical exponent. The presence of two length scales results 
in two types of dimensional crossovers. 



1. Quantum crossover 

A dimensional crossover in the quantum critical fluctu- 
ations occurs upon approaching the QCP either by lower- 
ing the temperature T or decreasing the tuning param- 
eter r; these crossovers are indicated by the horizontal 
and vertical shaded regions in Figs. [1] and [2] 

In the case of the 2d to 3d AF crossover, the positions 
of these crossover lines are given by T ~ A* and r ~ A^", 
where z = 2 and the correlation-length exponent has the 
mean-field value v = 1/2. For example, the specific heat 
coefficient at criticality, r = 0, changes its temperature 
dependence at the crossover temperature, T ~ A?, from 

7 ~ log(l/T) to 7 ~ const — y/T, as expected for critical 
thermodynamics of 2d and 3d AF fluctuations, respec- 
tively. 

The situation is more complicated for the 2d FM to 
3d AFM crossover, due to the mismatch of dynamical 
exponents in the two regimes. The low dimensionality, 
d = 2, of the FM spin fluctuations combined with a large 
dynamical exponent, z = 3, result in strong thermody- 
namic signatures of the 2d regime, which dominate over 
the 3d AFM fluctuations in an unexpectedly wide regime. 
For example, at criticality, r = 0, the 2d FM fluctuations 
yield singular quantum critical thermodynamics down to 
a temperature scale T ~ A*, where z = 3. Below this 
temperature scale, the thermal activation of these 2d fluc- 
tuations only yields Fermi-liquid behavior, but this is still 
much stronger than the contribution from the 3d part 
of the spin-fluctuation spectrum. Only at a much lower 
temperature scale, T ~ A*, the contributions from the 3d 
AFM fluctuations finally take over. Thus, we have here 
the peculiar situation that a 2d non-critical Fermi-liquid 
regime is sandwiched between the 2d and 3d quantum 
critical regions. 



2. Classical crossover 

In addition, there is a dimensional crossover associated 
with classical criticality. Upon approaching the classical 
phase transition line, T c (r), the correlation length £ in- 
creases. If £ reaches 1/A n , the classical fluctuations (i.e. 
those associated with zero Matsubara frequency) effec- 
tively change their dimensionality from 2d to 3d. This 
classical dimensional crossover occurs within the 2d quan- 
tum critical regime at the shaded line labeled T c i(r) in 
Figs. [1] and [2J and it also causes thermodynamic signa- 
tures. For example, the thermal expansion a(T) has a 
maximum at T c i for the 2d FM to 3d AFM crossover. At 
higher temperature T x , the crossover line T c ; enters the 
Ginzburg regime of the classical transition, and the clas- 
sical dimensional crossover becomes non-perturbative. 



3. Phase boundary and QCP location 

The phase boundary, T c (r), of the classical transition 
is linear in the distance to the QCP at elevated tempera- 
tures (with logarithmic corrections), but curves towards 
the QCP in the 3d regime. As a consequence, an extrapo- 
lation of the quasi-linear phase boundary in the 2d regime 
towards zero temperature yields an incorrect position for 
the QCP, see Fig. and we estimate the corresponding- 
error between the extrapolated and the true position of 
the QPT. 

We note that our treatment of thermodynamics is lim- 
ited to the non-ordered side of the phase transition and, 
in particular, breaks down upon entering the Ginzburg 
regime of the classical critical transition indicated by the 
thin line T G (r) in Figs. H and HI 

4- Observables 

As detailed below, we find that the thermal expan- 
sion and the compressibility are well suited to detect a 
dimensional crossover in the spin-fluctuation spectrum. 
Both possess pronounced signatures close to the expected 
crossovers as a function of temperature in the quantum 
critical regime, either a sharp drop or even a maximum. 
In contrast to this, the specific heat coefficient only shows 
a leveling off upon entering the 3d regime, which is harder 
to identify experimentally, see Figs. [5] and [5] below. 

B. Outline 

The body of the paper is organized as follows: In 
Sec. [n] we introduce the Landau-Ginzburg- Wilson field 
theory for magnetism near quantum criticality. We dis- 
cuss the anisotropic spin susceptibility and the associ- 
ated crossover in the Landau damping. We introduce our 
model for the dimensional crossover and give the result- 
ing formulae that determine the correlation length and 
other thermodynamic properties. Sees . IIIII and llVI are de- 
voted to a detailed discussion of the 2d - 3d AFM and 2d 
FM - 3d AFM crossovers, respectively. We shall derive 
phase diagrams and full crossover functions for thermo- 
dynamic quantities. Finally, in Sec.[V]we discuss existing 
experimental data vis-a-vis our theory results. We focus 
on the heavy- fermion metals CeCug-iAu^, YbRli2Si2, 
and CeCoIns which indeed display unconventional quan- 
tum criticality that may originate from quasi-2d spin fluc- 
tuations. A brief outlook concludes the paper. 

II. ORDER-PARAMETER FIELD THEORY FOR 
SPATIALLY ANISOTROPIC SPIN 
FLUCTUATIONS 

In order to analyze the dimensional crossover we will 
use the standard LGW critical theory of Hertz, Millis 
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FIG. 3: Schematics of a layered crystal structure with a weak 
hopping of spin-fluctuations between the planes, t' <C t. 



and Moriya for a (commensurate) itinerant paramagnet. 
The action of the Hertz- Millis- Moriya model reads^ 



S = l dT I ddr ^^o'HV,^)^ 



(1) 



where the real bosonic order-parameter field $ repre- 
sents commensurate spin fluctuations with a 3d ordering 
wavevector Q. We will generalize the field $ to have N 
components; the Heisenberg paramagnet corresponds to 
N = 3. The dynamics of the fluctuations are encoded 
in the propagator x^ 1 . Its momentum dependence will 
reflect the spatial anisotropy of the spin-fluctuation spec- 
trum. Its form will be motivated in the following. 



A. Bare susceptibility 

Starting from a model of interacting electrons on a 3d 
anisotropic lattice, a Hubbard-Stratonovich transforma- 
tion allows to introduce collective-mode variables repre- 
senting the spin-fluctuations. Their dispersion can be 
estimated, e.g., by RPA. Consider for simplicity a para- 
magnon dispersion of tight-binding type. On a 3d tetrag- 
onal lattice, a generic form is 

w(k) = t (2 — cos k x a — cos k y a) + 1' (1 — cos k z a) (2) 

where t and t' parameterize the hopping of the spin fluc- 
tuations within and perpendicular to the xy planes, see 
Fig. [31 the momentum k is measured relative to the or- 
dering wavevector Q, and a is a lattice constant. The 
quasi-2d character of the spin fluctuations is reflected in 
a small ratio between the hopping amplitudes, 



77 = *7*< 1. 



(3) 



Near the transition we can restrict ourselves to energies 
much smaller than the (large) in-plane bandwidth oc t; 
this is equivalent to a continuum approximation w.r.t. 
the in-plane lattice coordinates. The result is 



ui(k) = t$ I fcjf + ^ (1 - cos k z a) 



(4) 



where £o = a/V^ and ku is the in-plane momentum. 
(More generally, £o is a microscopic length scale of or- 
der 1/kp where kp is the electronic Fermi momentum.) 

The character of the spin fluctuations depends on 
whether their energy is larger or smaller than the ver- 
tical bandwidth t' . In the energy range f « w « (, the 
paramagnon excitation energy is mainly accounted for 
by the in-plane momentum fen , and in this sense the spin 
fluctuations are effectively two-dimensional, 



w(k) 



for t' < uj < t (2d). 



(5) 



The restriction on the energy oj implies that the in-plane 
momenta are confined to the window <C kn -C A, 
where A ~ t^ 1 is a momentum cutoff, and we defined for 
later convenience the crossover momentum 



A„ = ^A. 



(6) 



For such large energies, the paramagnons are ef- 

fectively dispersionless in the direction perpendicular to 
the planes, such that fluctuations with all vertical mo- 
menta k z are equivalent. On the other hand, for energies 
of order t' and smaller, the vertical momentum can be 
resolved, and the spin fluctuations have a 3d character. 
For to -C t' a full continuum approximation is justified, 



w(k) « <£o ( fcjj + nkl J , for (3d). 



(7) 



In this 3d regime, the momenta are now restricted to 
k z <C A and kn «A,. In both the energy ranges ui <C t' 
and f « u « (, the limiting forms of the dispersion, 
and (JTJ) , respectively, are sufficient for the description of 
thermodynamics. These forms are independent of micro- 
scopic details, the only requirement being the existence of 
a dimensional crossover scale oc t'. In contrast, the pre- 
cise properties of the crossover itself depend on details of 
the band structure, i.e., Eq. Q, and are non-universal. 
As we are less interested in these non-universal details, 
we approximate the dispersion by 



for A, ( < fc|| < A 



2d 



r\k\ for J k?, + r]k? z < A r) 3d 



(8) 



We will use the dispersion ([8|) for the propagator of the 
spin fluctuations 



(9) 



e 



for A,, < fc|| < A 2d 



T k ' u 1 k^+rjkl for Jkf t +rik 2 z < A n 3d 



Here, Sq represents the bare mass of the spin fluctua- 
tions, and the dynamics of the spin fluctuations is en- 
coded in the Landau damping term |w n |/rk which will 
be discussed in the next subsection. 
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Varying So drives the system through the QPT, which 
occurs at S = <5o jC r- We define the control parameter r 
of the QPT via 



r = S Q - So.c 



(10) 



such that r = at the QCP. To calculate observables 
like thermal expansion, we shall assume that the phase 
transition can be tuned by changing the pressure p, which 
is reflected in a pressure dependence of the bare mass 
So(p). 



B. Landau damping 

The dynamics of the paramagnetic spin fluctuations is 
controlled by damping due to particle-hole excitations in 
the metal; this Landau damping term |w„|/Fk dominates 
over the undamped dynamical term cx uj 2 in the propaga- 
tor at low energies. The function Tk has to be obtained 
from an expansion of the particle-hole bubble of the host 
metal, for small energies and momenta close to the order- 
ing wavevector. As the damping is different for ferro- and 
antiferromagnets, the different energy regimes discussed 
above need to be distinguished. Note that in the follow- 
ing we assume that the Fermi wavevector is sufficiently 
large for fermions to cause damping, 2k F > |Q|. 

Let us start with the dimensional crossover from 2d 
AFM to 3d AFM. In both asymptotic regimes, the 
wavevector Q connects hot lines on the Fermi surface, 
hence the known result for the antiferromagnet, Tk = 
const, applies. As we are not interested in details of the 
crossover, we shall employ 



T 2 d for A„ < feii < A, 



2d AFM 



Tk = 



T 3d for . k 2 + r\k 2 < A„, 3d AFM. 



(11) 



The dimensional crossover from 2d FM to 3d AFM is 
somewhat more complicated. The ordering wavevector 
is Q = (0,0, Q z ). In the 2d regime, i.e., for energies 
lu 3> t' , we assume that not only the spin fluctuation 
spectrum, but also the fermionic dispersion is 2d. With- 
out dispersion in the vertical direction, we obtain for Tk 
the standard result for the ferromagnet, Tk oc fey . On the 
other hand, in the 3d regime we again have Tk = const 
for the AFM. Hence, we approximate Tk by 



r 2 d Co^n for A„ < fen < A, 



2dFM 



r k = 



3d 



for 



kl + r]kl < A,„ 3d AFM 



(12) 



An explicit evaluation of the Landau damping for FM 
fluctuations in the presence of an anisotropic Fermi sur- 
face is given in App.JXJ The result given there, Eq. (|A6[) . 
reduces to the limiting form of (p~2|) in the 2d FM regime. 

In principle, the damping coefficients in the two differ- 
ent regimes, T 2 d and r 3( j, differ from each other. How- 
ever, as they only determine the overall energy scale, 



we will set them, for simplicity, equal from now on, 
r = r 2 d = r 3c ; (see also the discussion in Sec. lIIEI below'). 



C. Correlation length 

The quartic interaction uo among spin fluctuations in 
([T]) will modify its bare susceptibility For our model, 
this modification is captured by an effective correlation 
length £, 




(13) 



for A v < fc|| < A 2d 



|w„| + fey + rjkl for Jk^ + rjkl < A„ 3d 

where in the 2d regime the dynamical exponent is cither 
z = 2 for AFM or z = 3 for FM fluctuations. From 
now on, we will employ dimensionless units, i.e., we set 
effectively the length scale £o = 1 and the energy scale 
r = 1. In addition, we use a unit volume V = 1. 

For d + z > 4, the correlation length £ can be obtained 
from a self-consistent perturbation theory in uq'- 



r 2 = *<> 



N 



(> 



-u 



r^x(k,icj„). 



(14) 



Only in a regime of 2d antiferromagnetic fluctuations, 
this formula misses logarithmic corrections to the corre- 
lation length. In order to capture these, we will later 
have to apply the renormalization group (RG). 

Substituting the sum over Matsubara frequencies by 
an integral over the real axis, we obtain 



„ , . N + 2 r lw du 



u I — coth — 
6 J 7T 2T 



dkk 



Luk 2 



(15) 



A, 2n (r 2 + k 2 ) 2 + (Luk 2 - Z ) 2 

IT f A " dkk 2 UJ 



where A w is an additional energy cutoff, and we intro- 
duced the two-dimensional quartic coupling u = uqA/tt. 
The factor A/n originates from the dummy momentum 
integration over the z component in the 2d regime. The 
contribution in the second (third) line is attributed to the 
spin fluctuations of effectively two (three)-dimensional 
character. Correspondingly, we define a 2d (3d) regime 
in the phase diagram plane where the second (third) line 
in Eq. (fT5|) dominates the correlation length. We obtain 
the following criterion: 



T>Af- 1} or r>A^- 1} 



T < A 



2(z-l) 

n 



and r < A 



20-1) 



2d regime 
3d regime 



(16) 



where z = 2 for 2d AFM and z = 3 for 2d FM spin 
fluctuations. 
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D. Thermodynamics 

From the susceptibility (TIB"]) , we can obtain the free 
energy of the critical spin fluctuations 



NT 



^2 l0gX 1 (k,«Wn). 



(17) 



k,u; n 



It will be convenient to absorb a factor A/tt in the units 
of the free energy, Fir/ A — > F. Doing so, the free energy 
in our approximation takes the form 



F cr = — — [ — coth — 

7T 2T 



dkk to 
■ arctan 



7T 



A 



»7 JO 



■ arctan ■ 



2tt ^- 2 + fc 2 



(18) 



The integral in the first (second) line originates from the 
2d (3d) spin fluctuations. 

From the free energy, we can compute thermodynamic 
properties. We will consider the specific heat, thermal 
expansion, Griineisen parameter and the compressibility. 
The specific-heat coefficient 7 is defined as 



7 = 



d 2 F 
dT 2 ' 



(19) 



The thermal expansion a measures the change in volume 
as the temperature is changed, 



Vdf 



1 d 2 F 
V dpdT 



V~dp 



(20) 



Using a Maxwell equation, we have re-written the ther- 
mal expansion as a derivative of entropy with respect 
to pressure. In principle, all parameters of the model 
(TT|) might be pressure dependent. However, it has been 
argued 1 ^ that the most important contribution comes 
from the pressure dependence of the parameter multi- 
plying the most relevant operator in the model. Close to 
the quantum critical point, this is the control parameter 
r of the transition. Near a pressure-tuned quantum crit- 
ical point, we can expand the control parameter around 
the critical pressure p c , r w (p — p c )/po, where po is an 
a priori unknown pressure scale. In this case, there is 
a contribution to thermal expansion that measures the 
change of entropy upon variations of the control parame- 
ter r. Choosing dimensionless units, this contribution is 
given by 



dF cr 
dTdr 



(21) 



The critical Griineisen parameter r cr is the ratio of crit- 
ical thermal expansion and specific heat, 



r. 



Tic 



(22) 



The compressibility k measures the change in volume as 
the pressure is changed with temperature held fixed, 



]_dV 
V~dp 



1 d 2 F 

vW 2 



(23) 



The quantum critical fluctuations contribute to n an ad- 
ditive term, in the following denoted by K cr - The pressure 
dependence of the control parameter r results in the fol- 
lowing contribution to K cr : 



d 2 F cl 
dr 2 



(24) 



in dimensionless units. 



E. Universality? 

Before presenting actual results, it is worth asking how 
"universal" we can expect them to be. This question has 
various aspects: (i) Are the crossover functions (for a 
specific observable) universal in the sense that they do 
not depend on microscopic details - at least in a certain 
well-defined limit? (ii) Do all observables display the 
same crossover scale(s)? 

For aspect (i) the answer is that full universality does 
not exist, as the quantum phase transitions under con- 
sideration are at or above their upper critical dimen- 
sion. Therefore, even if the ultraviolet cut-off A and the 
crossover scale A^ are well separated, i.e. 77 1 , the 
bare values of A and the paramagnon interaction u will 
influence the crossover functions. In particular, the in- 
teraction u is at the origin of the classical dimensional 
crossover lines denoted as T c i in Figs. [HandO 

Moreover, as discussed at length in Sec. Ill Al the di- 
mensional crossover itself is determined by microscopic 
details, i.e. the precise crossover form of the bare suscep- 
tibility. At this point, possibly existing magnetic inter- 
layer frustration enters i 1 4 Similarly, the reference energy 
scales of the 2d and 3d regimes, T2d and T$d, which we 
have taken to be equal for simplicity, may differ (by a 
factor of order unity), which leads to a shift of the two 
asymptotic regimes on the temperature axis w.r.t. each 
other. 

For some of the observables listed above, additional 
care has to be taken regarding the dependence on exter- 
nal pressure. As the pressure dependence of all micro- 
scopic parameters is smooth, singular contributions to 
thermodynamics usually arise only through the pressure 
dependence of r in the vicinity of the QCP. However, in 
our model, the 2d regime possesses another relevant pa- 
rameter, namely the anisotropy r\. If the anisotropy also 
depends on pressure, it will give important additional 
contributions to thermal expansion and compressibility 
in the 2d regime. In dimensional units, these contribu- 
tions are represented by 



dF cr 
dTd^ 



2 F a 
drf 



(25) 
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They can be best estimated by considering the free 
energy F with a momentum dependence for the spin- 
fluctuations given by Eq. (j4|), instead of the simplified 
expression for F, (fT5|) . In the 2d regime, we obtain that 
the pressure dependence of the anisotropy yields contri- 
butions proportional to the ones deriving from a pressure- 
dependent control parameter, aj? r oc a cr and oc K cr . 
In the 3d regime, on the other hand, aj?r an d ^cr are 
suppressed by additional powers of momenta compared 
to a CT and k cx , respectively, and are therefore negligible. 
Having this in mind, we omit a further discussion of these 
terms in the following sections. 

Given all these caveats, our calculations to be pre- 
sented below are nevertheless valuable, for various rea- 
sons: (a) they illustrate the general structure of the 
phase diagram with the quantum and classical dimen- 
sional crossover lines, (b) they show the existence of the 
interesting intermediate regime where 2d quantum fluc- 
tuations coexist with 3d classical fluctuations, (c) they 
show which observables are especially sensitive to the di- 
mensional crossover and how "broad" or "narrow" the 
crossover signatures are expected to be, and (d) they 
show where logarithmic corrections can dominate power 
laws. 

Let us briefly comment on aspect (ii), namely whether 
all observables display the same crossover scale. For 
thermodynamics, there is a single relevant momentum 
crossover scale A^, at least in our simple model with- 
out inter-layer frustration. The single momentum scale 
however translates into various temperature scales, and 
depending on the thermodynamic quantity of interest dif- 
ferent crossover scales might be of importance. For the 2d 
- 3d AFM crossover, there are two relevant temperature 
scales: the quantum dimensional crossover temperature 
T ~ A 2 , and the classical crossover temperature T c i, de- 
termined by the paramagnon interaction u. Whereas the 
specific-heat coefficient is insensitive to the T c i scale, the 
thermal expansion and compressibility show signatures at 
both scales, and their dimensional crossover appears to 
be broad due to the existence of the intermediate region, 
A 2 < T < T c i. For the 2d FM - 3d AFM crossover, we 
find even three temperature scales arising from the single 
momentum scale A^: a quantum dimensional crossover 
temperature T ~ A*, a temperature scale T ~ A^ di- 
viding 2d non-critical from 2d quantum critical behav- 
ior, and a classical dimensional crossover scale T ~ T c \. 
These multiple crossover temperatures result in rich ther- 
modynamics and render the extraction of power laws es- 
pecially difficult. 

Transport properties - which are not subject of this 
paper - can be expected to display even more compli- 
cated crossover behavior. In addition to the momen- 
tum scale A^, there are characteristic length scales for 
transport scattering processes. For coupled chains, it has 
been argued 2 ^ that the interplay of those length and time 
scales can render transport fully three-dimensional even 
in a regime where the thermodynamic behavior is Id. 
Hence, thermodynamic and transport crossovers do not 



have to coincide in general. 



III. DIMENSIONAL CROSSOVER: 
ANTIFERROMAGNET TO 3D 
ANTIFERROMAGNET 



2D 



For materials consisting of antiferromagnetic planes, 
which itself are weakly coupled in the third direction, 
a scenario of a crossover from 2d antiferromagnetism to 
3d antiferromagnetism is plausible. As mentioned in the 
introduction, such a scenario might be realized in the 
heavy- fermion metal CeCus-zAu^. In the following, we 
present an analysis of this crossover within the framework 
of the LGW model whose applicability to CeCu6-a;Au x , 
however, has been questioned^^ 



A. Renormalization group 

The Hertz model JTJ for AFM fluctuations in 2d is at 
its upper critical dimension, d + z = d£ = 4. Conse- 
quently, there are important logarithmic corrections in 
perturbation theory that have to be summed, e.g., with 
the help of the renormalization group (RG). The RG 
equations for the running mass S(b) and quartic coupling 
u(b) in the 2d AFM regime readti*2i£ 



85 
d log b 

du 
d log b 



= 26 



N + 2 
12tt 2 



uS . 



N + 8 
12tt 2 



(26a) 
(26b) 



Here, we have introduced the RG scale b, and the RG 
flow corresponds to a reduction of the momentum-space 
cutoff, A — > A/6. The one-loop correction to S, i.e., the 
tadpole diagram, can be expanded in a power series in 6: 
The constant term is finite and non-universal (i.e. cutoff- 
dependent), whereas the linear term is universal. The 
latter linear term is written as second term in Eq. (|26aD ; 
the constant term will be absorbed in the initial value 
of 6. Hence, the flow starts at b = 1 with the initial 
conditions S(b = 1) = r^d (that differs from 5q by the 
non-universal Hartree shift) and u(b= 1) = u = uqA/tt, 
being the effective 2d quartic coupling. Note that the 
prefactors in (|2"6")l have been evaluated at T = 0, as their 
temperature dependence is subleading. 

Integrating the RG equation for the quartic coupling 
we obtain 



u{b) 



I2tt 2 



f 



N ■ 



8 log 



bA 
A 



(27) 



where we have introduced the momentum scale A, which 
depends on the bare quartic coupling constant u, 



Y2ts z 



A = AeW+§> 



(28) 
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Using the running quartic coupling, the solution for the 
control parameter is easily obtained. Substituting 5(b) = 
r(b)b 2 , we get 



-(b) = 



T2d 



log 



N + 2 
N + 8 



(29) 



As we will see later, the parameter rid differs from the 
control parameter, r = r^d, of the quantum phase tran- 
sition by corrections of order A 2 . 



B. Correlation length 

The logarithmic scale dependence of the coupling and 
the control parameter has to be in calculating the corre- 
lation length. Hence, Eq. (fT"5)) is replaced by 



r 2 = r(A/A*) 



N + 2 



■(A/A*) 



o 7T Ja v 2vr (^-2 + k 2y + , 



N + 2 iru* f A " dcu / >A " dkk 2 



6 A 



v Jo 





du) 


'o 


7T 


du> 




7T 


Jo 



2T 



(30) 
1) 



u) coth 



The abbreviations A* and u* are defined through 



A* = maxff 1 ^,} (31) 



and 



u* = u(A/A v ) 



12tt 2 
N + : 



log- 



(32) 



the latter being the quartic coupling at the crossover scale 
A,,. Eq. ([3H)l can be understood as follows: The first line 
accounts for the 2d quantum fluctuations at T — which 
have been re-summed in r(b) using the RG, see Eq. (f2"9"| . 
Finite-temperature corrections to this 2d result can be 
treated perturbatively 2 ^ and are in the second line. The 
third line, finally, is attributed to the 3d fluctuations. 

In the following, the limiting behavior of the correla- 
tion length is discussed in detail. 



The logarithmic dependence can be understood as an in- 
cipient correction to the mean-field value of the correla- 
tion length exponent umf = 1/2. 

Quantum critical regime, T 3> rid- The limiting be- 
havior of the remaining integral is given by the small lo 
limit of the integrand, such that 



r 



r(A/A* 



7T N ■ 



2 N + 8 



(34) 



Hence, in the quantum critical regime, the correlation 
length is determined by temperature, £~ 2 ~ T, up to 
logarithmic corrections. The logarithms arc cither cut- 
off by the correlation length itself or by the dimensional 
crossover scale A^. In the limit £ -2 ^> A 2 , the correlation 
length is asymptotically given by 



r 2 = r(4 



7T N 

2N~ 



2 log log 
— T 



A 



for T 2 > A 2 ,. 

(35) 



The log log dependence in the numerator can be at- 
tributed to the 2d classical fluctuations, i.e, the 2d- 
like Matsubara zero mode, within the quantum critical 
regime. [Depending on whether the first or the second 
term dominates in (|35p. one can further distinguish two 
sub-regimes within the regime denoted as the renormal- 
ized 2d classical regime in Fig. [TJ] When the correlation 



length exceeds the crossover scale, i.e. 
logarithms are now cut-off by A^, 



£_- 2 < A 2 the 



r(A/A„) + - 



ttTV + 2 1°S W 



2 N + 8' 



for r 2 <A, 2 



(36) 



The crossover at the scale £~ 2 ~ A 2 is associated 
with the advocated dimensional crossover for the classi- 
cal critical fluctuations, where the Matsubara zero mode 
changes its character from 2d to 3d. At criticality, 
f2d = 0, this occurs at a temperature 



Tcl\ r2d =o 



A 2 



log log £• 



(37) 



Correlation length in the 2d regime, T ;§> A 
rid > A 2 V 



In the 2d regime as defined in Eq. (fTB")) , we can neglect 
the last line in the expression (f3"0"|) as it yields only a small 
correction. We further distinguish two sub-regimes. 

Fermi-liquid regime, T -C rid- In the Fermi- liquid 
regime, the correlation length is determined by the so- 
lution of the RG equation for the tuning parameter 



r 2 = 



rid 



(log- 



N + 2 
N + 8 



(33) 



This classical dimensional crossover temperature, T c ;, is 
logarithmically enhanced as compared to the quantum 
dimensional crossover temperature, T ~ A 2 , see Fig. [TJ 

The classical dimensional crossover occuring at T c \ be- 
comes non-perturbative when the line T c i(r) enters the 
classical Ginzburg region. Generally, our perturbative 
treatment breaks down sufficiently close to the classical 
transition at the Ginzburg temperature Tg, where the 
classical Ginzburg parameter is of order one, Q ~ 0(1). 
For 3d classical fluctuations, the Ginzburg parameter is 
given by Q — U£, where the classical quartic coupling is 
U = ttu*T/A v ~ T/(A v log-£-). Upon increasing tem- 
perature, the crossover line T c i(r) approaches Tq(t) and 
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enters the Ginzburg region at a temperature T x where 
T G and T cl coincide, see Fig. [TJ 



T x ~ A 2 , log 



A 

A,, 



(38) 



2. Correlation length in the 3d regime, T <C and 
r S d < A 2 V 

In the 3d regime, we can instead neglect the second line 
in (|30|) . We again distinguish between two sub-regimes. 

Fermi-liquid regime, T <C r$d-' In the Fermi- liquid 
regime, we can set T = in ([3D)) in order to obtain the 
leading estimate for the correlation length, 



r3d 



T2d 



( l0g ^) 



A 7 + 2 
N l-s 



(39) 



where C is a non-universal constant that depends on the 
chosen cut-off structure. The control parameter r = r^d 
obtains a shift from the contribution of the parame- 
ter regime where the spin fluctuations have developed 
their three-dimensional character; the consequences of 
this shift are discussed in Sec. IIII B 31 The temperature 
correction neglected in (|39|) is of Fermi- liquid type ~ T 2 . 

Quantum critical regime, T 3> r 3 d- In the quantum 
critical regime, the temperature corrections dominate the 
correlation length. Evaluating the leading behavior we 
obtain 



r 



T3d 



7T\ 3 / 2 N- 

2) « 3/2) W- 



T 3/2 



(40) 



where we used the explicit expression for u*, Eq. 
The quantum critical regime can again be further subdi- 
vided into regimes where either of the above two terms 
dominates the correlation length. 



the non-perturbative character of the classical dimen- 
sional crossover. Neglecting the logarithmic corrections, 
the phase boundary varies linearly with temperature, 
T c oc rid, in the 2d regime. 

In the 3d regime, on the other hand, we use expression 
(|30|) and get 



T3d 



3/2 TV 

am 



t;' 



3/2 



AT + SA.log £ 



(42) 



The different behavior of the phase boundary in the 2d 
and 3d regimes, (|4"Tj) and (|42|) . implies that an extrapo- 
lation of the phase boundary from the 2d regime leads 
to an incorrect position of the quantum critical point. 
The 3d spin-fluctuations shift the quantum critical point 
slightly towards the ordered phase. The quantum critical 
point is not at rid = 0, but rather at rsd = 0, leading 
to a difference between the extrapolated and the actual 
position of the QCP, 



Ar = r 3d \ r2d=a 



A 2 

v 



log- 



(43) 



The position of the extrapolated QCP, Ar, see Fig. [4] is, 
however, located close to or even within the 3d pocket of 
the phase diagram. 



C. Thermodynamics 

We turn to a discussion of the thermodynamic quanti- 
ties specified in Sec. Ill Dl They can be obtained with the 
help of expression (TT5)) for the free energy, with the corre- 
lation length given by the self-consistent Eq. (|3H)) . A nu- 
merical solution for the specific heat, thermal expansion, 
and compressibility in the quantum critical regime is pre- 
sented in Fig. [5J A detailed discussion of their asymptotic 
behavior is given below. 



3. Phase boundary 

The position of the phase boundary, T c (r), in the (T, r) 
plane can be estimated by setting the correlation length 
to infinity. In the quantum 2d regime for temperatures 
lower than the crossing temperature T < T x , Eq. 
the classical dimensional crossover from 2d to 3d is per- 
turbative as it is located outside the classical Ginzburg 
region. Here, we can use expression (|3"o) for the correla- 
tion length to obtain 



TT N + 2 



rid 



2 N 



N + 2 
' N + S 



(41) 



For larger T > T x , the logarithmic corrections to the 
critical temperature T c will differ from Eq. (|4"Tj) due to 



1. Thermodynamics in the 2d regime, T>A, or rid S> A^ 

Within the 2d regime, the thermodynamics is domi- 
nated by the first line in the expression for the free energy 
Eq. (fTS]) . We again distinguish between two sub-regimes. 

Fermi-liquid regime, T -C r^d'- The evaluation of the 
specific heat and thermal expansion is straightforward 



7c 



N T 
12^' 



The resulting Griineisen parameter reads 

r =- l - 

2r 2dl0g^= 



(44) 



(45) 



The evaluation of the compressibility is more involved. 
Here, we have to take into account the effective momen- 



tum dependence of the correlation length, £ k 



r(A/fc), 



10 




extrapolated 
QCPs 



QCP 



control parameter r 

FIG. 4: Locating the QCP based on experimental data re- 
quires an extrapolation of the finite- T phase boundary to 
T = 0. Such an extrapolation procedure is prone to a sys- 
tematic error if a dimensional crossover at lower temperatures 
takes place. The extrapolation of the phase boundary (thick 
line) from the 2d regime towards T — yields estimates of the 
position of the QCP that are shifted towards the non-ordered 
phase by an amount Ar. The different thin lines and the as- 
sociated extrapolated QCPs illustrate that the extrapolation 
itself is ambiguous as the 2d phase boundary is not strictly 
linear in T but has logarithmic corrections. For the 2d - 3d 
AFM crossover, both the size of the 3d pocket and the shift 
Ar, Eq. (J43J), are of order A^. For the 2d FM - 3d AFM 
crossover of Sec. lIVI the size of the 3d pocket is instead given 
by A*, while the shift Ar oc A^, Eq. (|70|l . Consequently, the 
extrapolated QCP is outside the 3d pocket in this case. 



see Eq. (f2T)|) . in the expression for the free energy (fl~5|) 
that arises after RG improvement of perturbation the- 
ory. We obtain 



N f dk fdr(A/k)\ 



47T 2 ./a« & V 9r 2 d J 



(46) 



The asymptotic behavior originates from the lower limit 
of the momentum integral that has to be cutoff either by 
the crossover scale or the correlation length itself, A* = 
maxf^ 1 , A r) }. In the 2d regime we thus obtain for the 
asymptotic behavior; 23 ! 24 



N 



Air 2 1 



q N+2 

N+8 



loe 



A 



n JV + 2 



(47) 



Note that, for N < 4, Eq. ([47j) predicts a correction to 
the compressibility that diverges upon approaching the 
quantum critical point, T2d — > 0. This divergence is only 
cutoff upon entering the 3d regime. 

A diverging electronic compressibility has interesting 
consequences. In particular, a coupling of the electronic 
system to lattice degrees of freedom can render the cou- 
pled system unstable, resulting in a first-order transition 
driven by quantum critical AFM fluctuations. 25 Such a 
fluctuation-driven first-order transition may occur in the 
present case for a sufficiently large 2d regime. 




10 10 10" 10 

temperature T 



0.06 




10 



10 10 10" 

temperature T 



10 



FIG. 5: (color online) Crossover behavior of thermodynam- 
ics for r — 0, i.e., in the quantum critical regime of the 
2d AFM - 3d AFM crossover scenario. The vertical dashed 
line at higher T ~ T c i indicates the classical and the one at 
lower T ~ Aj| the quantum dimensional crossover tempera- 
ture. The insets show the same data on a double-logarithmic 
scale. The chosen parameters are A^ = A = 20, A,, = 0.01, 
U = 1 and N = 3. Panel (a): Specific-heat coefficient 7. 
Also shown are the asymptotic behaviors given by Eq. (I48p 
(red/dashed) for the high-temperature 2d regime and Eq. (|58[) 
(blue/dash-dot) for the low-temperature 3d regime. Panel 
(b): Thermal expansion a. The asymptotic behaviors are 
given by Eq. (|49)l (red/dashed) and by Eq. ((59} (blue/dash- 
dot). Panel (c): Compressibility n. The asymptotics are in 
Eq. d52j) (red/dashed) and Eq. JBD} (blue/dash-dot). Note the 
very weak increase with T in the high-temperature 2d regime. 
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Quantum critical regime, T r^d •' In the 2d quantum 
critical regime, the specific heat depends logarithmically 
on temperature, 



N A 



7c 

For the thermal expansion we obtain 
N d£- 2 (, T 



877 dr2d 



log 



A 2 



const 



(48) 



(49) 



This leads to slightly different asymptotic behavior 
within the 2d and 3d classical regime, see Fig. [TJ 



N + 2 
' N + 8 



log log jj- if T» 



N + 2 

N + 8 rp 

log w 



if 



T<T d 



(50) 



For T 3> T c i , we obtain for the Griineisen parameter the 
temperature dependence 



log log 



Tlo ^( lo 4) 



N + 2 
N + 8 



(51) 



The critical part of the compressibility is dominated 
by the Matsubara zero mode. Its form also differs in the 
classical 2d and 3d regime, 



dr 2 d 



N T 

N_ T 

16 A„5 



if t»t c/ 

T if T«T C/ 



(52) 



Using the expressions for the correlation length (|35[) and 
(|36p . we obtain the asymptotic behavior at criticality 

T2d = 0, 



— 9 A 



log 



log log 



vT 



A 



if T>T C , 
if T«T c; 



(53) 



The asymptotics of thermal expansion and compress- 
ibility in the logarithmically small intermediate regime 
A 2 < T < T c i are not shown in Fig. [5] 



2. Thermodynamics in the 3d regime, T< and 
r 3 d « A 2 V 

In the 3d regime, the critical contributions to thermo- 
dynamics are coming from the second line in Eq. (|18l) . 

Fermi-liquid regime, T -C r^d'- Here, the specific heat 
has the form 



N , A 

7 = "7T log -7- + 7cr, 

6 A„ 



7c 



12 A„ 



(54) 



The critical part, 7 cr , that depends on the control pa- 
rameter, r 3c i, leads here only to a small correction to 
the background contribution that originates from the 2d 
fluctuations. The thermal expansion reads 



Nit T 
24 A vy ^' 



(55) 



We omit here and in the following multiplicative factors 
that are powers of d£~ 2 /dr2d, that includes the logarith- 
mic normalization of the 2d control parameter, see (|39[) . 
In the Fermi-liquid regime the critical Griineisen ratio, 
r cr = a cx /(Tj CI ), reads 



1 1 

2^? 



(56) 



As explained in Ref.ll2L in the Fermi- liquid regime scaling 
predicts a universal critical Griineisen ratio in the sense 
that the proportionality factor in the relation r cr oc 1/r 
is just determined by critical exponents. The prefactor 
— 1/2 in Eq. (|56p is in agreement with this scaling pre- 
diction. For the compressibility we get 



N y/rjd~ 
Stt A„ 



(57) 



Quantum critical regime, T 3> r$d-' Again, the critical 
part of the specific heat is only sub-leading ([54")) now with 



(58) 



15C(5/2)ATVT 
7cr ~ v^32 A~~ 

The thermal expansion is given by 
a cr = N 



3C(3/2) VT 



2tt16 A ?) ' 



(59) 



The thermal expansion behaves as v T in the 3d regime. 
The sub-leading correction to (|59|) vanishes as T 3 / 4 , with 
a two-fold origin. First, there is a contribution due to 
the next-to-leading term in the expansion of the second 
derivative of the free energy d 2 F CI /dTd^ 2 that is of 
order £~ l /A v . Second, and more importantly, there is 
a T-dependent correction originating from the deriva- 
tive <9£ _2 /<9r 3 d that contributes to a a term of order 
T 3 ' 2 /£ _1 AH. Both contributions originate from the T 
dependence of the correlation length, induced by the 
bosonic interaction u£L The critical Griineisen parame- 
ter deriving from (|58|) and ([59]) obeys scaling, T CT ~ 1/T. 

The compressibility is dominated by the Matsubara 
zero mode, 



N T£ 
16 AT 



log - 



, T l/4 if 



/rldKi 



if 



r 3 d < 



r 3 d > 



r 3/2 



A„ log 



r 3/2 



A, log t; 

(60) 
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where we used the expression (|40j) for the correlation 
length. There exist two sub-regimes where the compress- 
ibility varies either linearly with T or, very close to the 
quantum critical point, as T 1 / 4 (Ref. \2w. Note that the 
critical contribution to the compressibility, K cr , changes 
sign upon crossing over from the Fermi liquid 15T|) to the 
quantum critical regime (|60|) . 



we obtain 



IV. DIMENSIONAL CROSSOVER: 2D 
FERROMAGNET TO 3D ANTIFERROMAGNET 

We now turn to the dimensional crossover from a 2d 
ferromagnet to a 3d antiferromagnet, describing the sit- 
uation of weakly AFM coupled ferromagnetic planes, 
where the 3d ordered state corresponds to so-called A- 
type antiferromagnetism. 

The key difference to the crossover in Sec. [Hi] is related 
to the form of Landau damping and, as a consequence, 
the different dynamical exponents z = 3 and z = 2 in 
the 2d and 3d regimes, respectively. The combination of 
low spatial dimensionality, d = 2, and large dynamical 
exponent, z = 3, results in strong thermodynamic con- 
tributions of the 2d FM spin fluctuations. Even below 
the temperature T < A 3 , thermally excited 2d fluctua- 
tions lead to a large non-critical Fermi liquid background 
that dominates over the critical 3d AFM fluctuations in 
a wide parameter range, see Fig. [2j This leads to a pecu- 
liar situation at criticality, r = 0: the 2d and 3d quantum 
critical regions are separated by a temperature regime 
A 4 < T < Af ( where non-critical Fermi liquid behavior 
prevails. 

Other differences between the present FM - AFM 
crossover and the AFM - AFM crossover of Sec. Ifffl 
are: For the FM - AFM crossover, the uniform sus- 
ceptibility changes from critical to non-critical behavior, 
see Sec. IIV CI Finally, the effective dimensionality is al- 
ways larger than the upper critical dimension, d + z > 4. 
Therefore, we can use directly Eqs. (|15[) and (Ti~g|) for 
the correlation length and free energy, respectively, with- 
out the need of a RG improvement. This simplifies the 
analysis considerably. 



A. Correlation length 

In the analysis of the asymptotic behavior of the cor- 
relation length, we again distinguish between the 2d and 
the 3d regime, where £ is dominated by the first and 
second integral in (fl"5|) . respectively. 



rid 



ir(N + 2) T 
-u 



144 



3/2 
r 2(/ 



(61) 



where r 2 d differs from the bare mass So by a constant 
cutoff-dependent shift. Here, the temperature-dependent 
part is of the Fermi-liquid form and subleading. 

Quantum critical regime, T > r 3 / 2 and T > A 3 : In 
the quantum critical regime, on the other hand, temper- 
ature dominates the correlation length 



N + 2 
24tt 



uTlog 



1 



(62) 



with A* = maxff^A,} as above, Eq. (J3TJ. The 2d 
quantum critical regime can be subdivided into three 
regimes depending on how £~ 2 compares with A 2 , and 
whether the first or the temperature-dependent second 
term dominates. At ~ A 2 , a dimensional crossover 
associated with the classical critical fluctuations takes 
place. At criticality, r w rid = 0, this crossover occurs at 
a temperature of order 



Tri 



u log ' 

Ji T7 



(63) 



At this classical dimensional crossover, the logarithmic 
increase of (|62|) with increasing correlation length £ is 
cutoff. This allows, in particular, a solution for the phase 
boundary, see below. However, at a higher temperature 



T 

-*- r 



A 2 



(64) 



T c i(r) enters the Ginzburg regime and the classical di- 
mensional crossover and, as a consequence, the expres- 
sion for the phase boundary becomes non-perturbative. 

Non-critical Fermi liquid regime, T <C A 3 and r 2 d "C 
A 2 ; The 2d critical behavior crosses over into a 2d non- 
critical Fermi liquid behavior below a temperature T ~ 
A 3 or, as a function of the control parameter at r ~ A 2 , 
see the shaded lines in Fig. [21 In this intermediate regime, 
the correlation length has the form, 



T2d 



V 



(65) 



where the numerical factor C\ is non-universal, i.e., de- 
pends on our modeling of the dimensional crossover. It 
turns out that this 2d Fermi-liquid background still dom- 
inates thermodynamics down to a temperature and con- 
trol parameter scale A 4 where finally the 3d critical be- 
havior takes over. 



1. Correlation length in the 2d regime, T ^> A 4 or 
r 2 d > A 4 

Fermi-liquid regime, T -C ri^ 2 and r 2d S> A 2 : Evalu- 
ating the correlation length in the 2d Fermi- liquid regime 



2. Correlation length in the 3d regime, T -C A^ and 
r 3 d < A 4 

Fermi-liquid regime, T <C r^d- The contribution from 
the 3d fluctuations shift the value for the control param- 
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ctcr 



r3d = r 2 d + C 2 uA r 



ir(N + 2) T 2 



144 



A 1 > 2 
IV V r 3d 



(66) 



where C 2 is a constant dependent on the employed cutoff 
structure. The temperature dependence is again of Fermi 
liquid type. Comparing the temperature corrections to 
the correlation length in the various Fermi liquid regimes, 
one obtains the two control parameter crossover scales, 
and A 2 that are shown in Fig.[2j 



Quantum critical regime, 
critical regime we get 



T r 3c i: In the quantum 



r 2 = r 3d + (N + 2) 



C(3/2) T 3 / 2 
u- 



24V2^ A r , 



(67) 



A comparison of this temperature dependence with the 
one of Eq. (|65j) yields the quantum dimensional crossover 
temperature, T ~ A^. 



3. Phase boundary 

In the limit of vanishing correlation length, Eq. 
yields for T c < T x , see Eq. (jM)) , the following expression 
for the phase boundary within the quantum critical 2d 
regime 



N ■ 



r 2 d = -- 



24tt 



-uT c log ; 



a 2 ' 



(68) 



It varies (up to logarithmic corrections) linearly with T. 
In the intermediate non-critical Fermi liquid regime, the 
phase boundary behaves as r 2 d ~ — uT c 2 /A 3 , before cross- 
ing over into the 3d regime. There, we use expression (|67p 
to obtain 



r 3 d 



-(N 



2)^ T ' 



3/2 



24\/2T A, 



(69) 



Similar to Sec. IIIIB 3) the extrapolation of the phase 
boundary within the 2d regime towards T = leads to an 
error in the estimate for position of the quantum critical 
point of 



Ar 



r 3 d 



T2d=0 



(70) 



B. Thermodynamics 

Thermodynamic quantities follow from the expression 
([18]) for the free energy. A numerical solution for the spe- 
cific heat, thermal expansion and compressibility in the 
quantum critical regime is shown in Fig. O In the follow- 
ing, an analysis of the asymptotic behavior is presented. 
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FIG. 6: (color online) Crossover behavior as in Fig. [5] but 
now for the 2d FM - 3d AFM case. The vertical dashed 
lines denote the three crossover scales T c i > > A^, see 
Fig. [2] T c i denotes the classical dimensional crossover where 
£ ~ 1/A^; T ~ separates the 2d critical regime from the 
2d non-critical Fermi-liquid regime; and below T ~ A* 3d 
quantum critical behavior sets in. The chosen parameters 
are A w = A = 20, A,, = 0.01, u = 1 and N = 3. Panel 
(a): Specific- heat coefficient 7, with the asymptotics given 
by Eq. (|74|1 (red/dashed) for the high-temperature 2d regime 
and by Eq. (158[) (blue/dash-dot) for the low-temperature 3d 
regime. Panel (b): Thermal expansion a, showing a max- 
imum at the classical dimensional crossover. The asymp- 
totics at high (red/dashed) and intermediate T (green/dash- 
dash-dot) are given by (|76|) . and at low T (blue/dash-dot) 
by Eq. (I59j> The inset shows that the asymptotic \fT behav- 



ior only sets in for T < A v 



10" 



Panel (c): Compress- 



ibility k, with the asymptotics in Eq. 1)78 p (red/dashed and 
green/dash-dash-dot) and Eq. (|60| (blue/dash-dot) showing 
the T 1 / 4 behavior that again only sets in for T < A 4 . 
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Thermodynamics in the 2d regime, T> A* or r^d A 



Fermi-liquid regime, T <C r^J d 2 and r 2 d 3* A 2 : In the 
2d Fermi-liquid regime, evaluating the leading behavior 
is straightforward. For the specific-heat coefficient and 
thermal expansion we obtain 



Ntt 1 



7c 



Ntt T 

'id 
12 



This gives a universal Griineisen parameter 

1 1 



r cr = 



2r 2d ' 



(71) 



(72) 



For the critical part of the compressibility we find 



Using Eq. ([62]) for the correlation length we obtain the 
asymptotic behavior for the compressibility at criticality, 
r « r 2rf = 0, 



-r-t-r. if T> Tci 



if T«T C( 



(79) 



Note that k depends singularly on the dangerously irrel- 
evant interaction u. 

The results ([76]) and ([79]) show that not only the 2d 
quantum critical regime above T c \ , but also the interme- 
diate 2d quantum 3d classical critical regime below T c i 
displays a well-defined universal asymptotic behavior. 

Non-critical Fermi liquid regime, T <C and r 2 d <C 
A 2 : In the intermediate non-critical regime, thermody- 
namic quantities display a Fermi-liquid form: 



(73) 



3/2 



and T > Ah In 



Quantum critical regime, T 3> r 2d 
the 2d quantum critical regime, the specific heat is given 
by 



(74) 



The thermal expansion depends logarithmically on tem- 
perature 



N T 2/3 

8^ l0g max{£- 2 ,A2} 



const . 



(75) 



with the correlation length given in Eq. ([6"2"]) . Using its 
explicit form, it becomes clear that the thermal expansion 
has a maximum at the dimensional crossover, £ _1 ~ A^, 
that occurs at a temperature T c ; and is associated with 
the classical Matsubara zero mode: 



log 



T 2/3 

log^p 



if T»T C/ 
const, if T<T C/ 



(76) 



In the quantum critical 2d regime, the Griineisen param- 
eter shows the asymptotic behavior 



1 



y2/3 : 



(77) 



where we omitted logarithmic corrections that depend on 
the effective dimensionality of the classical fluctuations, 
see Eq. ([76]). 

Similarly, the compressibility also exhibits an addi- 
tional dimensional crossover at T c i as it is determined 
by the classical fluctuations 




16 A,^- 1 



if T»T C , 
if T«T C/ 



1 r 



r T 



(80) 



with constants C3 and C4 being non-universal. The com- 
pressibility is determined by 3d classical fluctuations, 



_ N T 
Kcr ~ 16 A^T' 

where the correlation length £ has the form ([65 



(81) 



2. Thermodynamics in the 3d regime, T< and 
r 3d < A* 

The thermodynamics in the 3d regime resembles the 
one discussed in Sec. lIII Cl The background contribution 
for the specific heat coefficient in Eq. ([54"]) is, however, 
now determined by the 2d FM fluctuations, 7— 7 C 



1 

2A„ 



C. Uniform susceptibility 

A particular property of the 2d FM - 3d AFM 
crossover is that the susceptibility, measuring the re- 
sponse to a uniformly applied magnetic field, is critical 
within the 2d regime, whereas it becomes non-critical at 
the dimensional crossover. We consider the static limit 
and study its temperature dependence 



Xu (T) = x(T, a n = 0, fci, = 0, k z = Q z ) = 



1 



(82) 



As before, the wavevector k is measured relative to the 
(3d) ordering wavevector Q. £q 2 is a constant propor- 
tional to the hopping t' between the 2d planes, £g oc 
t'[l — cos(Q z a)]. From the discussion in Sec.[II]it follows 
(78) that £q 2 ~ A 2 . The temperature behavior is fully ac- 
counted for by the T dependence of the correlation length 
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FIG. 7: Behavior of the uniform susceptibility Xu f° r r — 0, 
i.e., in the quantum critical regime of the 2d FM - 3d AFM 
crossover scenario. The chosen parameters are £q 2 = 0.0001, 
A = 20, A„ = 0.01, 



A 



1 and N = 3. The vertical 



dashed lines represent the crossover scales Tq = £q ~ A v < 
A^ < A^. The dashed/red line is the 2d high-T asymptote 
(|83|) . and the dashed-dot/blue line is the 3d low-T asymptote 



£. In the quantum critical 2d regime, i.e., at high tem- 
perature, we can neglect the constant contribution £q 2 , 
and we obtain with the help of (|62]) for the asymptotic 
behavior 



Xu(T) 



Tlogi 



(83) 



This divergence is cut off below a scale Tq = £q 2 ~ A 2 , 
i.e. essentially at the classical dimensional crossover line. 
For T <C Tq, the constant contribution from the t' dis- 
persion dominates, Xu = const. In the 3d quantum crit- 
ical regime, the temperature dependence of the corre- 
lation length (|67|) leads to the T-dependent correction 
Xu =const +5xu(T), 



(84) 



in agreement with the results of Ref. 28j. The result of a 
numerically evaluation of the uniform susceptibility (|82p 
is shown in Fig. [7] The lower two crossover temperatures 
leave essentially no trace in Xu , as the regimes only differ 
in the small thermal correction to the large Xu(T = 0). 



V. APPLICATION TO HEAVY-FERMION 
METALS 

Our results are of potential relevance to layered nearly 
magnetic metals where indications for two-dimensional 
criticality have been found. In this section, we discuss 
the cases of the heavy-fermion materials CeCug-zAuz, 
YbRh 2 Si2, and CeCoIn 5 . While CeCoIn 5 possesses a 
layered lattice structure and should naturally display a 



dimensional crossover, for CeGvLQ- x Au x and YbRli2Si2 
there is only empirical evidence (discussed below) for 
quasi-2d critical behavior, with no obvious microscopic 
reason. We note that, in all three cases, ingredients be- 
yond the LGW spin fluctuation theory may be important 
for a full understanding of the critical behavior. 

We refrain from discussing strongly correlated 
transition-metal compounds with layered structure, such 
as high-temperature superconducting cuprates. In these 
materials, a standard LGW approach alone cannot be 
expected to capture the relevant physics due to the prox- 
imity to the half- filled Mott insulator. 



A. CeCus-xAuj 

CeCug is a paramagnetic heavy Fermi liquid, which 
can be driven into an antiferromagnetic metallic phase 
by Au substitution.— >2£ In CeCna^ x Au x , the quantum 
critical point is located at x c ~ 0.1. For x > x c , the 
AF order is known to be 3d, which implies non-vanishing 
magnetic couplings in all spatial directions. In contrast, 
neutron scattering in quantum-critical CeCus-zAu^ has 
provided direct evidence for a quasi-2d antiferromagnetic 
fluctuation spectrum. 6 This suggests that a dimensional 
crossover occurs near the quantum critical point, which, 
however, has not been experimentally identified to date. 
Let us therefore apply our results from Sec. IIIII in order 
to search for experimentally measurable consequences of 
a dimensional crossover from 2d AFM to 3d AFM. 

First, there is the location of the QCP itself. Ex- 
perimentally, the phase boundary appears to be linear, 
Tn oc (x — 0.1), in other words, the data points marking 
the finite-temperature phase transition can be linearly 
extrapolated to a putative QCP at x c = 0.1. From the 
consideration in Sec. Illll it is clear that the asymptotic be- 
havior of the phase boundary should be Tn oc (x — x c ) 2 ^ 3 , 
which suggests that the true x c > 0.1. Hence, samples 
with x = 0.1 may not be located at the quantum crit- 
ical composition. As a consequence, the system would 
be located outside the 3d pocket of the phase diagram 
in Fig. |4l this could then explain why no signatures of 
3d spin fluctuation were found in the neutron scattering 
experiment of Ref. 0. Unfortunately, concrete predic- 
tions are problematic due to the presence of logarithmic 
corrections in the d = z — 2 theory for the 2d AFM: In 
fact, in the 2d regime the phase boundary should fol- 
low Tn logT^ oc (x — x c ), Eq. (|4ip - this makes a linear 
extrapolation ambiguous. Experimentally, such logarith- 
mic corrections are difficult to extract. 

Second, the dimensional crossover should be manifest 
in thermodynamics. Whereas the crossover signatures in 
the specific heat are weak, the thermal expansion and 
the compressibility show a pronounced step-like behav- 
ior in the quantum critical regime of the 2d — 3d AFM 
crossover, see Fig. [5] The latter two are therefore more 
appropriate to detect a dimensional crossover in the crit- 
ical spin-fluctuation spectrum. 
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Thus, we propose to search for the dimensional 
crossover in CeCug-zAux by (i) detecting the change in 
behavior of the thermal expansion at criticality and (ii) 
by looking for deviation of the phase boundary from lin- 
ear behavior by employing pressure tuning of, e.g., an 
x = 0.2 sample. From the available data, the dimen- 
sional crossover temperature is likely below 100 mK. We 
also note that the presence of quenched chemical disor- 
der may modify the behavior near criticality at very low 
temperatures 



B. YbRh 2 Si 2 

The heavy-fermion material YbRli2Si2 displays a phase 
transition at 70 mK. 7 The low-temperature ordered 
phase is believed to be antiferromagnetic, although con- 
firming neutron scattering data are not available to date. 
The ordering temperature can be suppressed by applying 
a small field, resulting in a field-driven QCP at B = 0.06 
T in the ab plane and 0.66 T along the c axis. The or- 
dering temperature can also be suppressed by Ge doping: 
YbRh2(Sii- a; Ge 3 ;)2 with x = 5% seems to order at only 
20 mK in zero field. — The critical properties of YbRli2Si2 
are inconsistent with the predictions of LGW theory for 
3d AFM spin fluctuations. As for CeCu6_a;Au x , it has 
been speculated that the Kondo effect breaks down at 
quantum criticality - this idea received support from Hall 
effect measurements which indicate a pronounced change 
in the low-temperature Hall coefficient across the critical 
fields 

Remarkably, YbRli2Si2 appears to be almost 
/erromagnetic. This is particularly striking for 
YbRli2(Sii_ 2; Ge 2 ;)2 where the uniform susceptibil- 
ity follows Xu{T) cx T-°- 6 above 0.3 K. 8 In addition, 
the unexpected observation of an ESR signal below 
the Kondo temperature^ 3 , has been related to strong 
ferromagnetic correlations.— 

It is therefore worth discussing which properties of 
YbRli2Si2 appear consistent with the scenario of a 
crossover from 2d FM to 3d AFM, as would arise in a 
system of ferromagnetic layers with weak antiferromag- 
netic inter-layer coupling. 

Interestingly, thermodynamic measurements are par- 
tially consistent with 2d FM criticality, but below 0.3 K. 
The specific heat follows C(T)/T oc T~ ' 3 below 0.3 K, 31 
and the Gruneisen ratio diverges as T(T) oc T~ 0,7 be- 
low 0.6 K (Ref. [35h - these two exponents are close to 
the values —1/3 and —2/3 expected for 2d FM fluctu- 
ations. In addition, the temperature-field scaling ob- 
served in YbRh 2 (Sii_ a; Ge-^a with x — 5% (Ref. l36l) is in 
accord with ferromagnetic criticality, provided that one 
interprets (B — B c ) (where B c = 0.027 T is tiny) as the 
field conjugate to the order parameter. However, other 
observations in YbRli2Si2 appear inconsistent with this 
idea of 2d FM criticality: 37 for example, the fractional 
exponent observed in the T dependence of the uniform 
susceptibility Xu(T) cannot be easily explained with this 



scenario. 

In summary, the physics of YbRli2Si2 cannot be ex- 
plained in a straightforward manner in terms of near- 
critical 2d FM fluctuations (which turn to 3d antiferro- 
magnetism at lowest temperatures) alone. However, the 
experiments, showing critical ferromagnetic fluctuations 
which are cutoff only at very low temperatures, hint that 
a dimensional crossover of the type considered here may 
be important. Further investigations of YbRh.2Si2 sam- 
ples with Ir or Co doping 38 will shed more light on the 
role of the various crossover scales in this interesting ma- 
terial. 



C. CeCoIn 5 

The compounds CeMIns (M = Co,Rh,Ir, also dubbed 
115-compounds), unify a variety of fascinating heavy- 
fermion phenomena: CeCoIns and Celrlns are (likely 
unconventional) superconductors with T c = 2.3 K and 

0. 4 K, respectively, whereas CeRhIn5 is an antiferro- 
magnetic metal with Tn ~ 3.6 K. Transitions between 
these ordered phases may be tuned using pressure, chem- 
ical substitution, or magnetic field. In contrast to most 
other heavy-fermion materials, CeMIns is quasi two- 
dimensional, i.e., its lattice structure consists of weakly 
coupled layers. Consequently, a dimensional crossover 
scenario should naturally apply. 

A particularly interesting transition occurs in CeCoIns 
upon application of a magnetic field. 39 Superconductivity 
survives up to a critical field of H C 2 « 4.95 T. Normal- 
state properties near H C 2 are suggestive of quantum crit- 
ical behavior: both specific heat and resistivity display 
non- Fermi liquid temperature dependencies at H C 2, and 
the A coefficient of the resistivity diverges upon ap- 
proaching H C 2 from above. These features have been 
interpreted as signatures of an antiferromagnetic quan- 
tum critical point at (or close to) H C 2, with the ordered 
phase for H < H C 2 being suppressed by the onset of su- 
perconductivity. By applying hydrostatic pressure, the 
two phenomena - superconducting H C 2 transition and ap- 
parent magnetic quantum criticality - can be separated, 

1. e. H C 2 decreases faster than the quantum critical field, 
suggesting that the two phenomena are not related^ 

Recently, thermodynamic properties of CeCoIns near 
H C 2 have been studied in more detail^ At H = 5 T, 
signatures of a crossover at T* w 0.3 K between two dif- 
ferent singular behaviors were identified. This is best vis- 
ible in the thermal expansion, which was found to follow 
a(T)/T oc l/T in the temperature range 0.5 K < T < 
6 K, whereas a(T)/T is consistent with 1/y/T for 0.1 K 
< T < 0.3 K. The thermodynamic data below T* w 0.3K 
have been argued to be consistent with the predictions 
of the LGW theory for 3d AF spin fluctuations. (Note 
that the scale T* , below which 3d LGW behavior is seen, 
is shifted up to 1.4 K in doped CeCoIng_ x Sna..) For 
T > T* , the thermal expansion in CeCoIns seems con- 
sistent with 2d AF spin fluctuations of LGW type. How- 
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ever, the authors of Ref. |4l| have argued that the be- 
havior of the Griineisen parameter, being is somewhat 
reminiscent to that of YbRli2Si2, instead suggests non- 
LGW criticality (which nevertheless may arise from 2d 
critical magnetism). 

Taken together, the data indicate a crossover in crit- 
ical behavior at T* sa 0.3 K - this is also supported by 
resistivity measurements which show p{T) — p(0) oc T 3 / 2 
below T = 0.2 K at if = 5 T^ consistent with 3d LGW 
behavior. Although the nature of the critical behavior 
above T* is not fully understood - the presence of mul- 
tiple crossover scales^ 2 - complicates the analysis of the 
data - the interpretation of T* as a dimensional crossover 
scale is suggestive. The role of Sn doping in shifting 
this crossover scale is unclear at present; in a situation 
with geometric frustration of inter-plane magnetism one 
might speculate that disorder partially relieves frustra- 
tion. It would be worthwhile to investigate the magnetic 
excitations, e.g., of AF ordered variants of CeMIn 5 by 
neutron scattering, in order to determine the magnetic 
bandwidths in the directions parallel and perpendicular 
to the Celn planes. 



VI. CONCLUSION 
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We have studied the dimensional crossover of magnetic 
fluctuations in nearly quantum critical metals. Moti- 
vated by experiments on heavy-fermion systems, we have 
concentrated on the crossover from 2d FM or AFM fluc- 
tuations at elevated energies to 3d AFM fluctuations at 
low energies. Applying the standard Landau-Ginzburg- 
Wilson approach, we have obtained relevant crossover 
energy scales as well as crossover functions describing 
thermodynamic observables. 

The anisotropy in the spin-fluctuation spectrum leads 
to a dimensional crossover from 2d to 3d upon approach- 
ing criticality. We have found two types of dimensional 
crossover scales. Upon reducing temperature at critical- 
ity, there is a first dimensional crossover associated with 
the classical fluctuations and, at a lower temperature, a 
second crossover where the quantum fluctuations change 
their character from 2d to 3d, see Figs. [1] and [Jl In par- 
ticular, there is an extended intermediate temperature 
regime where 2d quantum fluctuations coexist with 3d 
classical fluctuations resulting in new power laws. For 
the 2d FM - 3d AFM crossover, there exists even a fur- 
ther sub-regime, where non-critical Fermi liquid behav- 
ior intervenes between the critical 2d and 3d regimes. We 
have found that the thermal expansion and the compress- 
ibility are well suited to detect a dimensional crossover 
of critical magnetic fluctuations. However, the existence 
of several dimensional crossover scales makes the exper- 
imental identification of asymptotic power laws particu- 
larly difficult. 



APPENDIX A: LINDHARD FUNCTION FOR A 
FERMI SURFACE WITH CYLINDRICAL 
SYMMETRY 

We evaluate the Landau damping of ferromagnetic 
fluctuations deriving from the Lindhard function for an 
anisotropic Fermi surface with cylindrical symmetry. The 
Lindhard function at T = is given by 

n(M"n) = -/^E . 1 . 1 • (Ai) 

J 2tt it - e q ml + iu n - e q+k 

In the following, we distinguish between a closed and an 
open anisotropic Fermi surface. 

1. Closed Fermi surface 

Here, we consider a closed anisotropic Fermi surface as, 
e.g., depicted in Fig. UK- Introducing the density of states 
and integrating over e, the dynamic part of the polariza- 
tion can be expressed as an integral over the (closed) 3d 
Fermi surface 

n dyn (fc,^„) = -iu n f — ^ — -, (A2) 

For an anisotropic Fermi surface, the density of states Vg 
and the Fermi velocity depend on the orientation of 
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the fermionic wavevector, q. We model the anisotropy of 
the Fermi surface with an in-plane Fermi velocity vp that 
differs from its z-component rjpvp. A large anisotropy 
with a quasi two-dimensional Fermi surface is obtained 
in the limit of small rjp. So we get 



1 /*2tt 

U dyn (k,iuj n ) = -— iuj n J d9 J < 



(A3) 



sin v v cos 



iuj n — vp(ht sin#cos</> + rjpk z cos (9) 



We further assumed that the density of states only de- 
pends on cos 9 with the azimuthal angle 9 of the Fermi 
momentum. In the limit of small w ra , this simplifies to 



IT dyn (fc, iu ri 



Vp 



{rj F k z 



(A4) 



where (v) is an angular average of the density of states 
that smoothly depends on momenta 



m-2 



dx 



Wo vT 



(A5) 



Generally, the damping of bosonic modes with momen- 
tum k is caused by particle-hole excitations close to the 
part of the Fermi surface that is tangential to k. This 
is reflected in the momentum dependence of {v). For 
the modes with vanishing in-plane momentum, k\\ = 0, 
only the part of the Fermi surface with azimuthal angle 
9 = ir/2 is involved in the damping processes, such that 
(i>) = u(0). On the other hand, the damping of modes 
with vanishing k z — can occur by exciting particle-hole 
pairs at any part of the Fermi surface and in this case 
(is) is a true average over the angle-dependent density of 
states. 

Thus we obtain for the damping function Tk, see 
Eq. ([9]), for an anisotropic Fermi surface 



r. 



2v F 



(A6) 



In the 2d regime, damping is dominated by the in-plane 
momentum, ~ |fcii|, and we obtain the expression ad- 
vertised in Eq. (|12l) . 



2. Open Fermi surface 



In the limit of a quasi two-dimensional Fermi liquid, 
the Fermi surface opens as depicted in Fig. [5Jd. If we 
neglect the warping of the Fermi-surface cylinder along 
the momentum q z direction, we end up with a Lindhard 
function of an effectively 2d Fermi system that is, in par- 
ticular, independent of the longitudinal momentum k z . 
Its dynamic part then has the form 



n dyn (fc||,iw„) = -v\\ 



v F k\\ 



E 



(A7) 



with the 2d density of states v\\ . The integration over the 
longitudinal momentum q z just yields a multiplicative 
factor. 

If we take the warping into account, we obtain instead 
the same limiting behavior as that of expression (|A4j) . 
In the limit rjpk z <C fcy, where rjp is again a small pa- 
rameter representing the strong anisoptropy of the open 
Fermi surface, the result is essentially unchanged from 
Eq. (|A7[) except that v» is replaced by an averaged den- 
sity of states. In the other limit of small in-plane mo- 
mentum rjpkz 3> fen, the warping in Fig.[8jD leads to a k z 
dependence of the dynamic part of the polarization re- 
sulting from particle-hole excitations now concentrated 
at the center, q z w 0, and at the edges, q z ~ ±n/a, of 
the Brillouin zone, where the Fermi surface is parallel to 
k z . 
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